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Two-loop contribution to the pion transition form factor vs. experimental 
data 

S. V. Mikhailov^*, N. G. Stcfanis 

'^Bogoliubov Laboratory of Theoretical Physics, JINR, 141980 Dubna, Russia 

We present predictions for the pion-photon transition form factor, derived with the help of light-cone sum rules 
and including the main part of the NNLO radiative corrections. We show that, when the Bakulev-Mikhailov- 
Stefanis (BMS) pion distribution amplitude is used, the obtained predictions agree well with the CELLO and the 
CLEG data. We found that no model distribution amplitude can reproduce the observed growth of the new 
BaBar data, though the BMS model complies with several BaBar data points. 



1. Form factor Y"' ^ in QCD, collinear 
factorization 

This transition form factor describes the pro- 
cess 7* (91)7* (92) — *■ T^^iv) and is given by the 
following matrix element {—q\ = > 0, — q| = 
g2>0) 

d4a:e-^^-^(7r«(p) | r{j^(z)j,(0)} | 0) = 

(1) 

For sufficiently large photon momenta Q^,q^ ^ 
m?p (where the hadron scale is set by the p-meson 
mass nip), all binding effects can be accumulated 
into a universal (twist-two) pion distribution am- 
plitude (DA), so that, on account of collinear fac- 
torization, one obtains the form factor as a con- 
volution: 



(low) factorization scale cannot be calculated 
from first principles and has to be modeled within 
some nonpcrturbativc approach, the amplitude 
T, describing the hard parton subproccsses, can 
be calculated in QCD perturbation theory. On 
the other hand, also the evolution of the pion DA 
with /ip is controlled by perturbatively calcula- 
ble evolution kernels and associated anomalous 
dimensions. Hence, in leading order (LO) of the 
strong coupling and at leading twist two one has 



pi 



dx- 



N 



f 



Q^x 



q^x 



^'•^\x)+0{Q-^) (3) 



with Nf ~ ^Jtt and x = 1 — x. The pion DA 



(2) 



is defined as 



0{Q- 



(2) 



|g(z)7^75C(z,0)g(0)|7r(P)) 



(2) 

While the pion DA ipir , which represents a pa- 
rameterization of the pion matrix element at the 
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(4) 



where C{z, 0) ^ V exp (ig ^p(r)c?T'^) is a path- 
ordered exponential to ensure gauge invariancc, 
and can be reconstructed from its moments, e.g., 
within the framework of QCD sum rules using 
either local [T] or nonlocal condensates |2l3j . 

The radiative corrections to the hard ampli- 
tudes T [Eq. ^] in next-to- leading order (NLO), 
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encapsulated in Ti, have been computed in [3|. 
More recently, the /3-part of the next-to-next-to- 
leading order (NNLO) amplitude T2, i.e., P^-Tp, 
was also calculated [5]. 

2. Light Cone Sum Rules for the process 

7*(Q2)7(q2c.O)^^o 

Experimentally, the transition form factor was 
measured by different Collaborations |6I7I8| when 
one photon is quasi real [q^ 0). This kinemat- 
ics requires the modification of the standard fac- 
torization formula Eq. ([2|) in order to take into 
account the long-distance interaction, i.e., the 
hadronic content, of the on-shell photon. To this 
end, Khodjamirian, [S], suggested a light-cone 
sum-rule (LCSR) approach, based on a disper- 
sion relation for F'^ in the variable q^: 



/>oo 

F^'"''^ {Q^,q^) = / 
Jo 



(5) 



The key clement of the LCSR is the spectral den- 
sity p(g2,s) = ^[F'^''^'''iQ^,-s)] for which 
one may employ the following ansatz [S]: p = 
pPi>(Q2, s)e{so -s)+ p^^[Q^,s)e{s - so), where 
the "physical" spectral density p^^ accumulates 
the hadronic content of the photon (below the ef- 
fective threshold sq) in terms of the form factors 
of vector mesons, viz., 

pP'XQ', s) - V2fvF^'^^Q^) ■ Sis - m\) , (6) 

where V stands for a p or an lo meson, 
while (F"^ contains the partonic part and is 
based on Eq. ([2|) via the relation (F'^ iCj)^ , s) = 
[(T ® ip^) {Q'^, -s)] . Using quark-hadron du- 
ality in the vector channel, it is possible to 
express the pion-photon transition form factor 
i^T*T"^(Q^O) in terms of p^'^: 

F^'^^^m = -r'^^^^^^^^^^^^^^ds 



1 r° Im(T(g)^^)(Q^-s 



)eK-^)A/^ds. (7) 



/o "V 

Here sq ~ 1.5 GeV^ and denotes the Borel 
parameter in the interval (0.5 — 0.9) GeV^. 

Partial results for p*^^^ at the NLO level have 
been presented in [TD], while the general solu- 
tion pi'^(g2,s) = IHI [(Ti(8^„)(Q^-s)] was 



recently obtained in [TT] : 

^J," =Cp|-3[l - «•(»)] + ^ 

-i„=(i)+2.(„)i„(i|)}.,:.„(x) 

n 

-Cf 2 ^ (G'„i + v{n) ■ hni) Mx) ■ 
1=0,2,... 



(8) 



Here {ipn} are the Gegenbauer harmonics which 
constitute the LO eigenfunctions of the Efremov- 
Radyushkin-Brodsky-Lepage (ERBL) evolution 
equation, with v{n),v"'{n) being the correspond- 
ing eigenvalues, whereas Gni, bni are calculable 
triangular matrices (see for details [TT]). 

Predictions for F'f'^'''{Q'^) at the NLO were 
given in [12] for various pion DAs, notably for 
the asymptotic one, the CZ model [1], and the 
BMS model [3]. We found that the result of the 
NLO processing of the CLEO data [7], the BMS 
bunch of DAs [3] , and the most recent lattice es- 
timates of the second moment of the pion DA are 
in good mutual agreement and inside the la error 

ellipse. The inclusion of the NNLO contribution 

(2) 

to the main partial spectral density p^ , propor- 
tional to /3o, was realized in [TT] taking recourse 
to the results of [5] • It turns out that it is negative 
and about -7%EI at small ~ 2 GeV^, decreas- 
ing rapidly to -2.5% at > 5 QeV^. The net 
result is a slight suppression of the prediction for 
the scaled form factor (see Fig. [T|). 

3. NNLO LCSR results vs. BaBar data 

Very recently, the BaBar Collaboration pub- 
lished new results on the pion-photon transition 
form factor that cover a wide range of momenta 
4 < < 40 GeV^ with high precision [8]. Sur- 
prisingly, their data exceed the asymptotic QCD 
prediction \/2/ir already at ~ 10 GeV^ and con- 
tinue to grow with up to the highest mea- 
sured momentum. This behavior would indicate 
that the 7*7 n process cannot be correctly 
described within the convolution scheme of QCD 
based on the coUinear factorization. We argued 



^takcn together with the effect of a more reahstic Breit- 
Wigner ansatz for the meson resonance in Eq. ^ 
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in [TT] that the inclusion of the NNLO radiative 
corrections cannot reconcile the BaBar data with 
perturbative QCD for any pion DA that vanishes 
at the endpoints a; = 0, 1 — see Fig. [T] Indeed, 
from Table [T] it becomes clear that also a wide 
pion DA. like the CZ model, cannot reproduce 
all BaBar data both in isolation or jointly with 
the CLEO data [7]. We conclude: 
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Figure 1. Q^F'''''^ (Q^) calculated with three different 
pion DAs: Asymptotic (lower solid line), BMS (shaded 
green strip), and CZ (upper solid red line). The BaBar 
data [8] are shown as diamonds with error bars. The 
CELLO [6] and the CLEO data [?] are also shown. The 
displayed theoretical results include the NNLO^ radiative 
corrections and the twist-four contributions. The hori- 
zontal dashed line marks the asymptotic QCD prediction 
V2U. 



(i) The combined effect of the negative NNLO^ 
radiative corrections and the Breit-Wigner ansatz 
for the meson resonances in Eq. ^ results in 
a moderate suppression of Q'^F'^ ^'^{Q^) in the 
range of momentum transfers 10-40 GeV^ [TT] . 

(ii) The hadronic content of the real photon is a 
twist-four contribution that is rapidly decreasing 
with increasing and can, therefore, not be the 
origin of the enhancement of the scaled form fac- 
tor measured by BaBar. 

(iii) Within the QCD convolution scheme, all 
pion DA models (cf . the in Table [T|) , which 
have a convergent projection onto the Gegen- 
bauer harmonics, and hence vanish at the end- 
points 0, 1, are in conflict with the BaBar data 
for g2F'''T"^(Q2) between 10 and 40 GeV^ (see 



Fig. [T]), because in this range these data violate 
the collinear factorization formula per se. 



Table 1 

for Asymptotic (Asy), BMS, and CZ DAs 



Pion 


BaBar 


BaBar 


BaBar only 


DA 


and 




10 data with 




CLEO 


all data 


Q2 > lOGeV^ 


Asy 


11.5 


19.2 


19.8 


BMS 


4.4 


7.8 


11.9 


CZ 


20.9 


36.0 


6.0 



4. Can the BaBar data be explained? 

There are no formal explanations of the BaBar 
data within the general framework of QCD at 
present — there are no reasons to adduce. How- 
ever, some theoretical scenarios have already been 
proffered |13|14|15] . We will discuss one class of 
such proposals based on the idea that the pion 
DA may be "practically flat" , hence violating the 
collinear factorization and entailing a (logarith- 
mic) growth of Q2^'^*7'r(Q2-) ^i^j^ q2 One has 

[H] (a = 0.53 GeV2; ^^{x) = A) 



3 Jo X 



1 — e 



dx . (9) 



Another option [T3] gives instead (m sa 0.65 GeV) 



V2 

3 



^ <^^(x, Q) 



dx 



X ■ 



(10) 



with ip^{x,^io) = U{N + (1 - N)6xx), {N w 
1.3). Eq. ([9]) can be compared with the available 
experimental data [6|7|8j by parameterizing them 
via the phenomenological fit (A w 0.9 GeV, b w 
-1.4) 



2V2/.^2 



A2 



A2 + Q2 



A^ 



A2 



as shown in Table [5J The main message is that 
one cannot fit the CELLO/CLEO data and the 
BaBar data with the same accuracy simultane- 
ously with both parameterizations. Now let us 



4 



S. V. Mikhailov and N. G. Stcfanis 



Tabic 2 

^ndf phenomcnological fit and the fit with a flat 

DA fike Eq. ^ — numbers in parentheses. The diagonal 
elements give the best-fit values of X^df ^^^^^ fix the cor- 
responding line parameters. 





CELLO& 
CLEO 


BaBar 


CELLO& 
CLEO 


0.48 (1.22) 


7.8 (15.8) 


BaBar 


10.8 (3.5) 


1.8 (1.8) 







push this point further and divide the BaBar data 
into two 'experiments' BaBarl and BaBar2, as 
indicated graphically in Fig. [2l One sees from 
Table [3] in terms of a Xndf criterion that the flat- 
DA scenario cannot describe both BaBar 'exper- 
iments' simultaneously with the same accuracy. 
One could interpret this outcome as an indica- 




Q'- [GeV^] 



Figure 2. SpUt BaBar data: ♦ - BaBarl, ■ - BaBar2. 

tion for an intrinsic inconsistency in the analysis 
of the BaBar data that deserves further attention. 

5. Conclusions 

In conclusion, we have presented a calcula- 
tion within the LCSR approach of the pion- 
photon transition form factor which includes the 
main NNLO radiative corrections and twist-four 
contributions. Our predictions are based on 
collinear factorization and agree with the data of 



Table 3 

^ndf ^^"^ flat-DA fit. Best-fit values on the diagonal. 





♦ BaBarl 


■ BaBar2 


BaBarl 


3.3 


0.33 


BaBar2 


3.5 


0.26 



the CELLO and the CLEO Collaborations but 
greatly disagree with the new high-Q^ BaBar 
data. Our analysis shows that all pion DAs, 
which vanish at the endpoints x = 0, 1, cannot 
reproduce the observed growth of the scaled form 
factor above 10 GeV^. On the other hand, flat 
pion DAs may describe this growth at high Q^, 
but at the expense that they fail to comply with 
the whole set of the BaBar data and also with 
those of the CELLO and the CLEO Collabora- 
tions. Future experiments may clarify this situa- 
tion. 
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